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Abstract 


We consider the Brownian motion of a cc pair produced in the very early satge of a 
quark-gluon plasma. The one-dimensional Langevin equation is solved formally to get 
purely mechanical properties at small and large times. Stochastically-averaged variances 
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1 Introduction 


Charmonium i.e., J/tp suppression [Q continues to be one of the most hotly debated sig¬ 
natures of the production of a quark-gluon plasma in ultrarelativistic heavy ion collisions. 
Both the initial formation and subsequent survival probabilities of the J/p) are affected 
by several factors viz. scattering with hard partons 0 , Debye colour screening |]I[ , drag 
and diffusion arising from Brownian motion [^, ^], evolution of the plasma via hydro- 
dynamic flow [Q, etc. Although the theoretical time scales for gluonic dissociation vs 
colour screened break-up are well known |P] the time scales for the swelling/ionization of 
J/ip caused by Brownian movement are not yet understood satisfactorily, and the aim 
of the present paper is to focus attention on this aspect. 

The classical 3-dimensional Fokker-Planck equation for the distribution function of 
a single charmed quark propagating in a plasma was first studied by Svetitsky |^. 
Assuming soft scattering with the partons he found the associated Boltzmann transport 
coefficients for drag and diffusion to be rather large but he ignored the force which binds 
the cc pair. Plotnik and Svetitsky ||^ extended this philosophy to the case of two-particle 
Fokker-Planck dynamics in the presence of a colour singlet/octet potential between the 
pair. However, since no attempt was made by them to actually solve the resulting 13 
variable partial differential equation, hence no simple estimate was given for the time 
scales of the cc pair. 

In the present work we adopt a different approach based on analyzing the one¬ 
dimensional Langevin equation for the stochastic trajectory of a cc pair initially bound 
by a screened Coulomb held. It is known since long ago that the Langevin theory 
provides a valid description of classical Brownian movement of a test particle acted upon 
simultaneously by a driving interaction, frictional force, thermal agitation, and random 
noise. In Sec.2 below we write formal solutions to the underlying equations of motion 
and obtain compact expressions for the purely mechanical observables at small and large 
times. Stochastic averaging of the observables is done in Sec.3 so as to deduce statistical 
properties {viz. means, variances, time scales, etc.) of the system. Sec.4 gives simple. 
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order-of-magnitude estimates of the relevant time scales and discusses the result vis-a-vis 
other mechanisms of J/ijj suppression. Finally, Sec.5 examines critically the validity of 
our main assumptions and also mentions several additional complications which would 
have to be incorporated in future applications of the theory. 

2 Purely Mechanical Observables 

2A. Assumptions & notations 

For a nonrelativistic cc pair propagating in a spatially homogeneous plasma the 
net external colour force due to the background is zero although the internal screened 
Coulomb potential (assumed to be colour singlet) survives. Working in the barycentric 
frame and adopting a one-dimensional view we can describe the motion of an effective 
single particle by dehning 

t = time, m = reduced mass, x = position, u = x = velocity, 

V = V{x) = binding potential, / = f{x) = —dV/dx = binding force, 

7 = coefficient of friction (per unit mass) assumed constant, 

T = ambient temperature of medium in energy units, 
h = h{t) = random force taken as a white Gaussian noise, 

C = 2m'-jT = strength of the noise autocorrelation function, 

D = = T/my = Einstein diffusion coefficient (2-1) 

The space derivative /' and total time derivative / of the force are written as 

/' = df/dx, f = df/dt = uf (2.2) 

Since the pair potential V depends upon |a;| the symmetry relations 

V(-x) = V(x), f(-x) = -f(x) (2.3) 

hold as shown schematically in Fig.l 
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Figure 1: Schematic shapes of cc binding potential V and force /. 

Before proceeding further a few important comments are in order. One-dimensional 
stochastic models P, || have been found very useful in the past because they are math¬ 
ematically simple and can also simulate purely radial motion in three dimensions. The 
choice C = 2m'yT guarantees that the test particle’s distribution at asymptotic time 
would become Maxwell-Boltzmann at background temperature T in accordance with 
the fluctuation-dissipation theorem |]T^ . Due to the same reason the single diffusion pa¬ 
rameter D = T/m'j gets hxed in terms of the temperature and the damping coefficient. 
All these remarks are, however, subject to alterations as will be pointed out later in Sec. 
5. 


2B. Langevin equation & the velocity 

The basic ordinary differential equation to be considered is 

ii + ■ju = {f + h)/m, x = u (2.4) 
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subject to the initial conditions 


u{t = 0) = Mo, x{t = 0) = xo (2.5) 

For a free particle and a harmonic oscillator, Eq.(2.4) was solved explicitly by Chan¬ 
drasekhar 1^ but the present case is more difficult because / is a nonlinear function of 
X. Using the integrating factor we obtain a formal solution for the velocity as 

u = +V + W (2.6) 


Here the contributions arising from free Rayleigh motion, the driving force, and the 
random noise are respectively given by 


u 


fr 


UqC ' , 


-it ft 

'■— / 
m Jo 


p-it ft 

w = - / dtie^^^hi 


(2.7) 


m JO 

with ti being an integration time and /i = / |ti, hi = h |ij. 

At small times the driving force can be approximated by a hrst-order Taylor expan¬ 
sion 


fi - fo + uofoh-{ -, 7ti < 1 


( 2 . 8 ) 


where the suffix 0 refers to the instant when the J/'ip was produced and the dots represent 
nonleading terms. Substituting into Eqs.(2.6), (2.7) we get the initial behaviour of the 
velocity and its square as 


uo + Jot-\ - \-w , t <C 7 ^ 

(2.9) 

{mq -|- 2uoJot -|- ■ ■ ■} -|- 


2{'Uo + Jot + ■ ■ ■}'zu -|- 

(2.10) 


with 


Jo = fo/m--fUo 


( 2 . 11 ) 
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Since the piece w containing noise may fluctuate rapidly with time its Taylor expansion 
is not attempted. Next, at large times the piece v in Eq.(2.7) can be integrated by parts 
once to yield the Rayleigh estimate 

nRi//m7 + ---, yt 3> 1 (2-12) 

But the asymptotic value of / is essentially zero both for a bound system (where the 
particle tends towards a point of stable equilibrium) as well as unbound one (where the 
particle tends to fly away). Therefore, we arrive at the leading asymptotic behaviours 

M ~ tc + ■ • •, ~ + ■ • •, t 3> 7“^ (2.13) 

2C. Analysis of Langevin trajectory 

Integrating Eqs.(2.6), (2.7) with respect to t we obtain the position, i.e., the relative 
separation between the cc pair 

X = + y + z (2.14) 


Here the contributions arising from free Rayleigh motion, the driving force, and the 
random noise are respectively read-off from 


fr 

X = Xo + 

7 (1 - -d. 


1 



y = — 1 

1 dtiKuJi, 


my J 

0 


1 

ft 

(2.15) 

z = - j 

I dtiKu^hi 

my J 

0 

with Ktti being a useful kernel defined by 

Ktt, = 1 - 

g-7(t-tl) 

(2.16) 


At small times the Taylor expansion (2.8) of the force can be used to deduce the following 
initial behaviour of the trajectory and its square 

X ^ XQ + Uot+'^t^+ + ■ ■ ■ + z, t < 7“^ (2-17) 

2 D 

X K 

x^ ~ {xq + 2xQUQt + (mq -|- XoJQ)t^ + {uqJq H- - — )t^ -!-■■■} 

+2{xo + Uot + -\ - }z + z‘^ (2.18) 

2 0 
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where 


j. fo^o 7/0 ,2 /o 1 n^ 

Ko = - h7Wo (2.19) 

m m 

fr 

At large times the sum x + y tends to the quantity 

IZ-n 1 

AToo = Xq -\ -1-/ dtifi, yt 3> 1 (2.20) 

7 my Jo 

whose value, however, is not known apriori. For a bound system may coincide with 
a point of stable equilibrium in the field /. However, if the particle becomes unbound 

r r 

then Xoo may look like Xq + w^/y with being the hnal velocity of free Rayleigh 
motion. Thus we obtain at the asymptotic behaviour 

X K. Xoo + 2:, 

x^ = + 2X^z + t > 7“^ (2.21) 

where no asymptotic expansion is attempted for the fluctuating term 2 :. 


2 D. Treatment of Langevin energy 

Finally we turn to the mechanical energy B = mu^/2-hV. Remembering the dash-dot 
notation specihed by Eq.(2.2) the rate of change of B becomes 

B = muii + V'u 

= {h — m'yu)u (2.22) 


whose formal solution is 


B 


Bn 


dti (jiiUi 


nxyul 


(2.23) 


At small times the initial behaviour (2.9 - 2.10) of the velocity can be inserted into 
Eq.(2.23) to yield 


B^ Bo 


[Uo + w)h — myluo 
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At large times if the system remains bound then E tends to a value below the ionization 
thresholds. However, if the system does ionize then the velocity u tends to w while the 
short-range potential V approaches zero. Hence, for disintegrated pair 


E ^ mtc^/2 f 3> 7 ^ 


(2.25) 


Eqs.(2.7 - 2.24) describe the main, purely mechanical, properties of interest to us; many 
of these expression may be regarded as new for general shape of the driving force /. 

3 Stochastic Averaging 

3 A. Statistical input 

In analogy with Chadrasekhar’s work ^ stochastic means will now be taken with 
respect to the initial velocity uq, initial position xq, and the Gaussian distributional of 
the noise h. For a genuine bound state obeying the virial theorem the input expectation 
values read 


(uo) = 0, (mq) = (^o) = 0 

{xD = A^o, {xqUq) = 0, {xofo)/m = -{uD 

{h,h) = C5ih-t2) 


(3.1) 


where A„o and A^-q are the velocity spread and position spread, respectively, in the 


barycentric frame of the J/ifj produced at f = 0. At general t the fluctuating velocity 
piece w and fluctuating position piece z have the properties 


(w) = 0, (z) = 0, 




(3.2) 


The informations (3.1) and (3.2) will be utilized below. 




3 B. Initial behaviour of averaged observables 

Let us return to Eqs.(2.9), (2.10), (2.17), (2.18), (2.24) describing the velocity, po¬ 
sition, and energy at small times. Clearly (Jq) = 0 and {Kq) = 0 since /o is an odd 
function of Xq in view of the assumption (2.3). Therefore, using the inputs (3.1 - 3.2) 
we deduce 


{u) = 0, Al = (n^) ^ + . . . 

(x) = 0, = {x^) ^ A^o + Ixt^ + --- 

{E) ^ {Eq) + g^t +■ ■ ■ , (3.3) 

where A„ is the velocity spread, A^, is the position spread, the dots stand for nonleading 
terms, and the extra coefficients written are 

9u = (C'-27m^Ay/m^ = 27(T/m-A^o) 

h = (C-f 27m(a;o/o))/3m^ = 5(„/3 
{Eq) = mA^o/2 + (Co), g^=mgu/2 (3.4) 

Eqs.(3.3 - 3.4) have fonr important physical conseqnences : 

(i) The linear time-dependence of the velocity variance A^ is controlled by the coeffi¬ 
cient gu- Clearly A^ increases with t if > 0 i.e. if T > mA^Q which is a nonequilibrium 
situation. However, A^ remains constant in time if g^ = 0 i.e. if T = mA^Q which is an 
equilibrinm sitnation. 

(ii) The cnbic time-dependence of the position variance A^ is governed by the pa¬ 
rameter Jj.. Evidently A^ increases with t if > 0 and the increment becomes 
comparable to the initial valne A^q at a time Tx satisfying IxT^ = A^q. In other words. 
Brownian movement can cause the bound state to swell snbstantially within a time scale 

Tx={AlJlxy^" < 7 -' if lx>0 (3.5) 

(hi) The linear time-dependence of average energy (E) is controlled by the constant 
g^^ = mgu/^. Obvionsly, (E) increases with t if g^>f) i.e. if T > mA^Q and it becomes 
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zero at a time satisfying = —{Eq). In other words, Langevin dynamics can 

cause cc to ionize after a time span 

= -{EQ)lg^ if > 0 (3.6) 

(iv) At this stage a couple of remarks must be added on the situation where inequalities 
on the coefficients get reversed, i.e. 

4 < 0 or fi-jj < 0 (3.7) 

The possibility Ja, < 0 physically implies that the contraction caused by the force term 
i^ofo) dominates over the expansion caused by noise C in Eq.(3.4) so that swelling of 
the bound state is ruled out. Next, the possibility < 0 in Eq.(3.3) implies that the 
mean energy (E) becomes deeper than (Eq) which again forbids the break-up of the 
classical bound state. Of course, the results contained in Eqs.(3.3 - 3.7) are new and 
original for a general shape of the binding force f{x). 


3C. Asymptotic behaviour of averaged observables 

Let us assume that J/ip has disintegrated under the above-mentioned conditions 
and examine Eqs.(2.13), (2.21), (2.25) describing mechanical properties at large times. 
Employing the statistical inputs (3.1 - 3.2) one finds 



~ -I- ■ ■ 

~ ‘lEt -\- 

<1 


~ T/m + 

{E) 

~ T/2 + --- , 

t 3> 7“^ 


where the contribution to arising from the variate Xoo has been omitted. Of course, 
Eq.(3.8) coincides with the well-known treatment of free-particle Brownian motion. Let 
us now apply numerically the results of the present section to the Langevin dynamics of 
J/ip produced in a quark-gluon plasma. 
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4 Numerical Work and Discussion 


4 A. Choice of parameters 

We use h = c = 1 units. The attractive potential and force between the cc are taken 

as g g 


V 

f 


4as 

3 I X 


exp 




V sign (a:) 


(4.1) 


where as is the squared QCD coupling constant, b the Debye screening length, and 
b~^ = Hd the screening mass. Various parameters of interest have typical values [§ 
given by the following two sets : 


Set I 


m = 0.75 GeV, Og = 0.4 
/i£) = b~^ = 0.2 GeV, 6 = 1 fm 

T = 0.2 GeV, 7 = 0.081 fm"^ (4.2) 

and Set II 

m = 0.75 GeV, a* = 0.6 

Hd = b~^ = 0.2 GeV, b = 1 fm 

T = 0.2 GeV, 7 = 0.157 fm"^ (4.3) 


Note that 7 has been identihed with the drag coefficient A appearing in Fig.2 of Ref. [Q 
and Fig.4.1 of Ref. |3 after replacing an incorrect numerical factor of 1024 by 512. 


4 B. Initial semiclassical properties 

Strictly speaking, the variances A^q and A^q should be obtained from the exact, 
3-dimensional, quantum mechanical, s-state Schrodinger wave function 'i/jq of the cc pair 
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at t = 0. However, for our simple, phenomenological purpose it will suffice to invoke the 
uncertainty principle for writing |mo| ~ h/m\xQ\. Then the semiclassical energy reads 

h;o = l/(2m|a;o|') + H(|a;o|) (4.4) 


Minimization with respect to |a;o| is achieved by setting dEQ/d\xQ\ = 0. This yields the 
condition 

Amas\xQ\ ( |a;o|\ ,, 

1 =- — -(6+IxoDexp I ——j (4.5) 

which can be solved numerically to get the size |a;o| in terms of b. Thereby we can 
estimate 


A, 


2 2 
^ Ur ^ 


uO 




Hn 


m‘^\xo\ 

\xq\^ as function of b 
' -b \ 1 


b + \xo\J m\xo\^ 


(Eo) - jul + Vo 


' b- |xo| \ 1 

,b+ |a;o|/ 2m\xo\‘ 


(4.6) 


4 C. Time scales 

It is now straight forward to evaluate the coefficients gu, Ix = and = mgu/^, 
dehned by Eq.(3.4). There are three time scales viz. 7“^, Tx = (A^q/Jj,) ' , and 
= —{Eo)/gj^ of interest in the present problem. All our numerical results are sum¬ 
marized in Table 1 corresponding to the two parameter sets (4.2 - 4.3). 


4 D. Discussion 

From a physical viewpoint the time scales 7“^, and represent respectively the 
frictional relaxation, positional swelling, and approach to ionization. A glance at Table 
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1 reveals that in the case of Set I (characterized by a weaker coupling constant o;* =0.4)) 
both Tj; and are positive and less than 7“^. Hence Brownian movement can cause a 
genuine break-up of the cc bound state in accordance with Eqs.(3.5, 3.6) above. 

On the other hand, in the case of Set II (characterized by a stronger coupling constant 
as =0.6), Tj; is imaginary and are negative i.e. unphysical. Hence random force plus 
diffusion cannot cause the cc bound state to dissociate in accordance with Eq.(3.7) above. 

Before ending we must remark that, in the context of J/ijj suppression, Langevin 
dynamics seems to be almost as important as other mechanisms (such as gluonic disso¬ 
ciation and Debye screening) invoked to explain the RHIC and LHC data. Evaluation 
of the charmonium survival probability [2,6] as a function of the transverse momentum 
reveals that typical time scales corresponding to the Debye and/or gluonic mechanisms 
are 5 - 10 fm/c. These numbers are quite comparable to the Langevin times Tx and 
of Table 1 (Set I) inspite of the differences in the input parameters T and iid- 


5 Additional Complications 


We now examine critically the validity of several oversimplihcations done above and 
also point out some additional complications likely to arise in future applications of the 
theory : 

(a) Bound state in 3-dimensions : 

One may argue that the simple, 1-dimensional, uncertainty principle based treatment 
of Eqs. (4.4, 4.5) will break down for the real charmonium which is a bound state in 3- 
dimensions. To answer this, we replace |a;o| by ro (which is the absolute distance between 
the cc pair) and appeal to the semiclassical, circular, Bohr orbits picture analogous to the 
familiar hydrogen atom problem. The ground state orbit in a Yukawa force / dehned by 
Eq. (4.1) has principal quantum number n = 1, orbital angular momentum Lq = h = 1, 
and centrifugal force condition 


mul Lq 

ro mrg 


(5.1) 


13 



which is entirely equivalent to the transcendental equation (4.5). It follows that our 
earlier estimate (4.6) of the bound state energy remains correct even for the real char- 
monium. 

(b) 3-dimensional random walk : 

Next, it is worth asking whether the main results of Secs. 2, 3 will get drastically 
altered if the random walk occurs in actual 3 dimensions. To answer this, we look at 
the vector Langevin equation 

(111 -* -* 

— + ^u={f + h)/m (5.2) 

where formal solutions, Taylor expansions, and stochastic averaging may be done on the 
same lines as Eqs. (2.6 - 3.2). Modifications appropriate to 3 dimensional configuration 
space can be readily done at every step. For example, at t = 0, we would have (/o) = 0, 
by parity argument while {fo.fo)/m = —{uq^) by the virial theorem. The crucial point 
to be noted is that the velocity variance A| will increase with t linearly and the position 
variance will do so cubically at small times, i.e., the essence of our leading behaviour 
(3.3) would remain intact. 

(c) Choice of initial conditions : 

One may claim that the time t = 0 should be set at the instant when the cc pair was 
created in the plasma by a hard partonic collision having divergent trajectories. In other 
words, the Langevin dynamics ought to have been applied even to the “pre-resonance 
formation stage” where some energetic pairs would lose their excess energy by random 
walk to form a bound cluster, not necessarily an s- wave ground state. This view, though 
very correct and ambitious, has three practical difficulties. First, pioneers like Xu et al 
and Karsch [2] have not adopted this view. Second, the initial values of 

(xo), {uo), (xqUo) (5.3) 

are not known immediately after the hard partonic reaction. Third, the final Brownian 
variances based on Fq. (5.3) will contain a large number of undetermined coefficients. 
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The initial conditions at t = 0 imposed in the present work correspond to a ^^fully - 
formed cc discrete bound state". This view, though modest, has three practical advan¬ 
tages. First, some pioneers of gluonic dissociation [2] have taken the initial state to be 
a standard cc resonance like fj', x etc. Second, the stochastic inputs (xq) = (mq) = 
(xqUq) = 0 are precisely known. Third, the hnal variances in Eq. (3.3) involve only one 
known effective coefficient Qu¬ 
id) Use of barycentric frame ; 

One may raise the criticism that there is no freedom to go to the cc barycentric system 
because the plasma - which is the source of random noise - provides a hxed frame of 
reference. For a plasma at overall rest this criticism is readily met by remembering that 
the noise h{t) being a function of time is Galilean invariant. Indeed, in a general frame 
of reference, the pair Hamiltonian H 12 reads 

Hi2 = P?/2mi -h pll2m2 + Uqgp + V{xi- X 2 ) - Xihi - 0:2/12 (5.4) 

where Uqgp is a constant potential generated in a spatially homogeneous plasma. Ob¬ 
viously, a separation between the centre of mass coordinate {xi + a^ 2)/2 and relative 
coordinate {xi — X2)/‘2. can be effected in Eq. (5.4). 

There is, however, an important word of caution here. In reality, the plasma evolves 
rapidly with the time by virtue of longitudinal/transverse expansion. The above-mentioned 
passage to the cc centre of mass frame is justihed in Bjorken’s boost-invariant hydrody¬ 
namics [6] if the cc pair moves either longitudinally or has small transverse momentum 
Pt- The procedure, however, may not be justihable for large pp pairs. This is because 
the pair distribution should relax to the equilibrium form in the plasma rest frame which 
would look different in other frames. 

(e) Miscellaneous rehnements : 

There are a few other subtle points to which attention will have to be paid in future. 
Since the cc pair may be in a colour singlet or octet state [7] a coupling between these 
channels may occur in the equation of motion. Next, the possibility of having different 
diffusion constants along the longitudinal and transverse directions should be allowed 
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so that the asymptotic equilibrium distribution of the charmed quark acquires a Tsallis 
shape [11] instead of the Boltzmann form. Finally, since the initial state of the char- 
monium is necessarily quantum-mechanical, a path-integral based density matrix may 
be formulated by taking hints from single-particle [12] or multiparticle [13] quantum 
stochastic dynamics. 
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Table 1: Numerical results on J/ijj suppression due to Brownian movement in the 
barycentric frame. For notations see the text. 


Set 

ko| 



(E,) 

7 ^ 

Tx 



(fm) 


(/™^) 

(GeV) 

(fm/c) 

(fm/c) 

(fm/c) 

I 

0.553 

0.226 

0.306 

-0.0245 

12.279 

5.175 

9.8447 

II 

0.346 

0.579 

0.119 

-0.107 

6.37 

(-3.64)^/^ 

-2.91 
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